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ABSTRACT 

Observations by the Goddard-University of Maryland plasma experiment 
on Explorer 34 have been used to study the fluid properties of the solar wind. 
By combining the plasma data with magnetic field data from Explorer 34 a 
study of discontinuity surfaces has shown that (a) the magnetic field vectors 
are parallel to some discontinuity surfaces as required for a tangential dis- 
continuity, (b) there are large bulk-speed discontinuities which probably indi- 
cate gliding motions at tangential discontinuities and are consistent with Sen’s 
stability conditions, (c) there is indirect evidence for the Kelvin-Helmholtz 
instability at some discontinuities in the solar wind. Several propagating 
hydromagnetic shocks observed satisfy the Rankine-Hugoniot conditions for an 
isotropic, single-component fluid. A study of the causes of sudden commence- 
ments and sudden impulses showed that ssc’s are nearly always caused by 
hydromagnetic shocks, but sudden impulses can be caused by shocks, hydro- 
magnetic tangential discontinuities or small dense regions in the solar wind. 
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On a much larger scale, it was found that the bulk speed and temperature of 
the solar wind are related by the equation ft = .036V - 5.54. The Sturrock- 
Hartle model gives results which are consistent with this relation, but applies 
only at very quiet times. The V-T relation can be described by an extended 
2-fluid model which postulates a heat source extending out to~20R . There 
are some regions in which the temperature is anomalously high, and they 
usually occur at positive bulk speed gradients and are probably due to turbu- 
lent heating produced in isolated patches by colliding streams; however, such is 
not the dominant heating process of the solar wind. Heating seldom occurs 
at negative gradients, suggesting that the Kelvin -Helmholtz instability may 
not be a dominant large-scale process. A study of the helium in the solar 
wind showed that it moves with the same speed as the protons, its density 
relative to the proton density is variable and often correlated with solar activ- 
ity, and its temperature may be 2 to 4 times greater than the proton tempera- 
ture suggesting heating at some point by the dissipation of hydromagnetic 
waves. The flow behind the shock of 11 August 1967 was studied using simul- 
taneous plasma data from Mariner 5 and Explorer 34. There was a high 
degree of correlation despite a separation of 1.6 x 10 7 km, indicating a quasi- 
stationary flow regime. The results suggest that the shock was driven by a 
broad, but not spherically symmetrical, high-speed plasma stream. 
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HYDROMAGNETIC OBSERVATIONS IN THE SOLAR WIND 


I. INTRODUCTION 

Experiments are conducted in the solar wind to determine the velocity 
distribution function of the plasma, which is a function of time, space, and the 
, species of particles studied. As has been pointed out before by Vasyliunas (1969), 

the length scale of the instruments used is so much less than any of the characteris- 
tic lengths of the plasma, that these instruments are just devices for the analysis 

V 

and measurement of charged particle beams, and do not exploit directly any of the 
cooperative properties of the plasma. What is done in practice is to determine 
the numbers of particles ~ dv in various differential regions of velocity space 
by counting them, and then either to fit to an assumed function, or use these 
values of differential density to approximate the unknown function. There are thus 
two problems, which have been fully discussed by Vasyliunas; 

1) The conversion of observed particle counts to values of ~ dv, this is 
the unfolding problem, which requires a good knowledge of the 
characteristics of the instrument, and 

2) The approximation to the distribution function 

v If we choose to restrict ourselves, or are restricted by lack of knowledge, 

to the fluid approximation, then we characterise the state of the plasma by a few 
* parameters which describe the distribution function and are analogous to 

statistics. In this approximation, the fluid parameters, density, velocity, and 
temperature are the statistics which describe its distribution function. When 
these quantities are supplemented with magnetic field measurements , one has 
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a complete hydromagnetic description of the plasma. These can be determined 
by fitting the differential densities to a known form of distribution function or by 
using the method of moments (Ogilvie, Burlaga and Richardson, 1967). The 
hydromagnetic approach is the appropriate one by which to study the situation in 
which there are scalar pressures for ions and electrons, large lengths, disturb- 
ances with rather low frequencies and very large electrical conductivity. 

It has been shown by Hundhausen et al. (1968) that the mean coulomb 
collision time in the solar wind becomes of the order of the expansion scale time 
at .1 AU so that collisions between charged particles cannot maintain an isotropic, 
maxwellian distribution function beyond 0.1 AU, and similarly collisions of ions 
with neutral hydrogen which may be present are also unimportant. The result is 
that the temperature parallel to the magnetic field tends to become greater than 
the perpendicular temperature, Tj_ as the plasma moves away from the sun. 
However, observations show that T I1 /T 1 ~ 2 at 1 AU (Hundhausen et al., 1967) 
where the magnetic field usually makes a 45° angle with the earth sun line, so 
an experiment which measures the temperature, T, along the earth-sun line gives 
a temperature which is between T„ and and usually differs from these by 
<50%. Since temperatures vary from < 10 4o K to ~ 10 6 °K, a factor of 2 is not 
important in many cases , and the fluid can be considered to be isotropic with 
temperature T to zeroth approximation. Presumably, this near isotropy is 
maintained in the absence of collisions by instabilities (Kennel and Scarf, 1968). 

The Goddard-University of Maryland plasma experiment which formed part 
of the experiment complement of the satellite Explorer 34 had an acceptance 
angle of 2.5° in azimuth and +9° in elevation, which was parallel to the satellite 
spin axis and perpendicular to the ecliptic plane. The time resolution is 3 
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minutes. Each spin, corresponding to observation at one energy per unit charge, 
was divided up into 16 equal angular intervals of 22 l A. degrees width. This angular 
resolution is useful in the study of the flow of plasma around the magnetosphere, 
but does not allow us to detect the thermal anisotropy of the plasma. In this way, 
Burlaga and Ogilvie (1968) examined the change in bulk speed across the earth's 
bow shock and the flow direction behind the shock along the flanks of the magneto- 
sheath, and found the results to support the fluid model. 

The fluid parameters were obtained by taking moments of the distribution 
function as briefly described by Burlaga and Ogilvie (1968), and the temperature 
is measured along the earth-sun line. 

The subjects to be discussed are as follows: 

II. Hydromagnetic Discontinuities in the Solar Wind 

III. Interaction of Discontinuities with the Earth 

IV . Solar Wind Heating 

V. Helium in the Solar Wind 

VI. The Flow Behind Shocks 

This choice of topics, and the discussion which follows, emphasizes our own 
work and interests but no attempt has been made to provide a comprehensive 
review of these subjects. 


3 


II. DISCONTINUITIES IN THE SOLAR WIND 


Two types of hydro magnetic discontinuities have been identified in the solar 
wind - tangential discontinuities and shocks. Burlaga (1968) showed that there 
are simultaneous discontinuities in the plasma and magnetic field parameters in 
the solar wind. He found several signatures which are consistent with those for 
tangential discontinuities but not for shocks. To prove that such a discontinuity 
is indeed tangential, one must show that 3 necessary conditions are satisfied: 

(1) it does not propagate relative to the solar wind, (2) the total pressure does 
not change across the surface, and (3) the magnetic fields B and B 1 on either 
side of the discontinuity surface are perpendicular to the surface normal, ft. 
Condition 1) was discussed by Fairfield (1968) and Burlaga (1968) showed a 
cascade of discontinuities which satisfied condition 2. Explorer 34 plasma data 
were used together with magnetic field data from Explorers 33, 34, and 35 to show 
that some discontinuities satisfy condition 3 (Burlaga and Ness, 1969). 

\ an ordinary fluid, relative motions along the surface of a tangential 
discontinuity give rise to the Kelvin-Helmholtz (K-H) instability and thus cannot 
persist. In a hydromagnetic fluid, however, relatively large motions are allowed 
because of the stabilizing effect of the magnetic field. Several very large dis- 
continuities in the bulk speed (>60 km/sec in < 3 min) were found in the 
Explorer 34 data and were attribut ed to relative motions along the surface of a 
tangential discontinuity (Burlaga, 1969a) Each of these discontinuities was 
associated with a directional discontinuity in the magnetic field, and the flow 
direction changed across at least 2 of the glide planes. The change in the 
magnetic field direction, w, tended to be near 90° and was never less than 46°; 
this is somewhat surprising, since Burlaga (1969b) showed that directional 


discontinuities tend to have small w. It was found that the large velocity dis- 
continuities with w near 90° were consistent with the conditions for stability with 
respect to the Kelvin-Helmholtz mode (Sen, 1963, 1964), but that similar dis- 
continuities with w near zero or 180° would be unstable. Thus, the unusual w 
distribution could indicate that large velocity discontinuities with small w dis- 
integrated by the K-H instability. 

Hydromagnetic shocks have been studied by several groups, but only the 
work of Sonett et al. (1966) and Ogilvie and Burlaga (1969) is based on simul- 
taneous magnetic field and plasma data. Our work concerns the Rankine-Hugoniot 
conditions for 6 propagating shocks. It was found that density jump and in most 
cases the temperature jump is given within experimental errors by the R-H 
conditions for an isotropic, single fluid plasma. The values of shock velocities, 
density ratios and Mach numbers So obtained indicate that at 1 AU the typical 
interplanetary shock is not strong. This agrees with the predictions of numerical 
calculations carried out by Hundhausen and Gentry (1969), and with the average 
observed transit time for shocks which can be associated with flares. 

Knowledge of the direction of the shock normal is necessary in order to 
apply the Rankine-Hugoniot conditions. Experience shows that this quantity is 
hard to determine from observations made by a single satellite, because the 
change in magnetic field direction is often small, or not accurately known. Our 
results include the description of shocks whose normal directions were found by 
the interpretation of magnetic field observations from more than one satellite. 

III. INTERACTION OF DISCONTINUITIES WITH THE EARTH 

A geomagnetic sudden commencement or sudden impulse can in principle 
be caused by a shock in the interplanetary medium whose front intersects the 
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magnetosphere, or by a discontinuity convected past the magnetosphere with the 
bulk speed of the plasma. In a frame at rest in the plasma, pr assure is balanced 
across such a discontinuity by definition; since it persists. The earth sees dynamic 
pressure changes at the passage of such a surface. 

Gold (1955) suggested that ssc’s were caused by shock waves propagating 
through the interplanetary medium from the sun. Subsequently, Sonett et al. 

(1966) reported direct interplanetary observations of a shock-like discontinuity 
which was moving from the sun and was associated with a ssc that was observed 
by 51 ground stations. Other observations of propagating events thought to be 
shocks, for which both plasma and magnetic field information was available, have 
been described by Ogilvie and Burlaga (1969). Nishida (1964) suggested that ssc’s 
could also be produced by a non-shock mode, presumably a hydrom&gnetic wave 
or a tangential discontinuity. 

Nishida (1964) also suggested that the non-shock mode discontinuity must be the 
cause of si, and Sonett and Colburn (1963) proposed that the si" is generally due 
to a reverse shock. Gosling et al. (1967) presented interplanetary observations 
which showed a discontinuous decrease in density and a gradual increase in 
temperature at the time of a world-wide, discontinuous decrease in the earth's 
field, thus establishing that an si could be caused by a non-shock mode dis- 
continuity. Ogilvie et al. (1968) examined simultaneous interplanetary plasma 
and magnetic field data associated with a similar si” and observed a discontinuous 
decrease in density, and discontinuous increase in magnetic field intensity and 
no appreciable change in temperature, showing that the si" was caused by a 
hydromagnetic discontinuity whose signature was that of a tangential discontinuity. 
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We have also shown that positive impusles si + are sometimes caused by 
hydro magnetic discontinuities. Gosling et al. (1967) reported an observation of 
a discontinuous decrease in temperature at the time of an ssc, showing that the 
event was not caused by a forward shock, but not ruling out a reverse shock. 

It is clear from these observations that both si and ssc can be produced in 
a variety of ways . 

F 

Taylor (1968) examined the causes of 36 ssc events using interplanetary 
magnetic field data from Explorer 28. He concluded that 8 of these events were 
due to tangential discontinuities, and that 26 were possible shocks, which caused 
the ’’larger" ssc events. Nishida (1964) suggested that the rise time of the 
impulse is small (£2 min.) for the events caused by shocks and large (>2 min.) 
for events caused by thicker, non-shock mode disturbances which propagate 
slowly or not at all. 

It seemed of interest to carry out a similar study using both interplanetary 
plasma and magnetic field data to see if it is possible to predict, using only 
ground observations, which type of interplanetary structure is responsible for a 
given si or ssc. 

A theory of the interaction of the solar wind with the earth shows that the 
change AH in the horizontal component of the earth’s magnetic field should be 
proportional to the change in momentum flux. Siscoe et al. (1968), from a study 
of 13 si + events, showed that from this theory one can calculate the change in 
momentum flux using A H and an empirical constant of proportionality which 
possibly varies with time. Ogilvie et al. (1968) showed a similar result for both 
si and ssc events. Clearly, A Ii alone is not sufficient to distinguish between a 
shock and a tangential discontinuity as the cause of a given event. We have 
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combined magnetic observations by Ness and Fairfield and our plasma observa- 
tions for 19 ssc’s and si’s during the period June to December 1967, in order to 
determine whether one particular type of hydromagnetic structure is associated 
with a particular type of event, and whether one can estimate the thickness of the r 

causative discontinuity from the rise time of an event as suggested by Nishida. 

As this work is described elsewhere (Burlaga and Ogiivie, 1969a), the 
procedure will be outlined. In order to classify the events as ssc or si, all 
events reported by 10 or more stations during this period were examined. A 
parameter A was defined using the number of stations calling a given event an 
ssc (#ssc) and the number calling it an si (#si) 

#(ssc) - #(si) 

A = 

#( ssc) + #(si) 

Thus if all the stations are unanimous in calling an event an ssc, A = +1, while 
if they all classify it as an si, A = -1. Table I shows details of the events and 
values of A. One can see from the table that there are basically two classes of 
event, but in many cases there was no concensus among observatories as to 
which a given event belonged. 

Examination of the data showed that all the events in Table I were associated 
with abrupt changes in the state of the solar wind near the earth. It has already 
been established by Ogiivie et al. (1968) and by Ogiivie and Burlaga (1969) that 
8 of the events were caused by hydromagnetic shocks; these events are indicated 
by the asterisks in Table I. Of the remaining eleven, four resembled shocks 
(July 25, Oct. 28, Nov. 3, Dec. 29) but were weak or ambiguous, five were 
probably tangential discontinuities and two were caused by small dense 
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regions in the solar wind, of spatial dimension ~.005AU, which produced 
corresponding 'pulses’ in the magnetogram. 

These interpretations are shown in the last column of Table I, where we 
see that the eight most prominent events (those reported by >32 stations) were 
all caused by hydromagnetic shocks. It is significant that not all of these events 
were classified as sudden commencements - for Aug. 29, A = -.19 and for 
Nov. 29, A = 0.1. If we define a sudden commencement by the criterion 
A > 0.75, then all sc's were caused by shocks, or conversely that A > 0.75 
indicates that the earth intersected a shock. 

The sudden impulses, defined by A < -.75, were caused by interplanetary 
structures, but not of a single type. Among the rather small selection of events, 
there appear to be two in which a sudden increase in the H component is followed 
within minutes by a sudden decrease to the prepulse value. 

Events in Table I with -.75 < A < .75 were not classified unambiguously as 
sudden impulses or sudden commencements. Five of these were caused by 
shocks, and three were caused by tangential discontinuities. Two of the three 
tangential discontinuities were associated with negative A and 3 of the 5 shocks 
were associated with positve A. There is thus a tendency for shocks to be 
identified as sudden commencements and tangential discontinuities as sudden 
impulses, but the relation is not good enough for predictions. 

When we examine the rise-time t given in Table I for the magnetic 
disturbances associated with these events, we see that the rise times for the 
seven shock events with A > 0.75, are not systematically shorter than those of 
events caused by tangential discontinuities. 
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In general, although the thickness of a discontinuity in the interplanetary 
medium is not related to the rise time of the corresponding event, and one 
cannot associate the rise time of a sudden commencement with the velocity of 
the causative shock, the disturbances which are observed in the geomagnetic 
field have their origin in characteristic disturbances in the hydromagnetic 
fluid-shocks and tangential discontinuities. 

IV. HEATING IN THE SOLAR WIND 

We now consider the temperature of the solar wind, and the heating 
mechanism by which it is produced. Figure 1 shows a distribution of the square 
root of hourly average values of T plotted against bulk velocity V . Although a 
good deal of scatter is evident, a functional relation is suggested. When this 
data is replotted by computing T from three hour averages for bulk speed 
intervals 250 km/sec 1 < V < 300 km/sec" 1 , etc, in order to take out short 
period fluctuations in T and V, we see the result in Figure 2. Here the present 
observations are shown as open dots, and the variability is shown by the error 
bar on the uppermost point. Various theoretical predictions are shown on this 
diagram (Parker's original isothermal model, and the results of various single 
fluid models) as well as certain other experimental results. The point marked 
W is from Wolfe's observations (Wolfe and McKibben, 1968) on Dec. 15, 1965. 
The average values of T x /i and V for the flight of Mariner n in 1962 are given by 
the point N (Neugebauer and Snyder, 1966). The point H represents observations 
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by Hundhausen referring to quiet times during 1962-1967. It is clear from the 
way these points lie on the line that it represents a valid relation between the 
average velocity and average temperature of the solar wind, in the sense that 
the solar wind properties are represented by a point on this line whose position 
varies with solar activity, but the slope of the line remains constant. 

The functional form of the relationship, 

T = (0.036 ±.003) V- (5.54 ±1.5) (1) 

is empirical (Burlaga and Ogilvie, 1969b). The relation V 2 = a + bT, which is 
suggested by the Bernoulli equation, gives a somewhat worse fit. 

It is widely believed that the predictions of the 2-fluid model, of Sturrock and 
Hartle (1966) are inconsistent with the observations, giving temperatures which 
are too low; this has led to the search for interplanetary heating mechanisms. 
However, Figure 2 shows that the predictions of Hartle and Sturrock (1968) are 
actually consistent with the extrapolated observations. Their 2-fluid model is 
simply a model for the very quiet solar wind, as the authors themselves have 
emphasized, and should not be applied to non-quiet conditions. Although the 
2-fluid model does give lower than average temperatures, it also give!* lower 
than average speeds. Since k (T e ± T. )«mV 2 , the low speed is a more serious 
problem than the low temperature. 

Parker (1963) showed that, if heat were added to the plasma such that its 

t 

temperature remained constant out to tens of solar radii with no heating beyond, 
high wind velocities could be obtained. Using Parker’s model, Burlaga and 
Ogilvie calculated T and V at 1 AU, assuming T = 10 6o K out to a radius R* . 

The results are shown in Table II. Such a crude calculation gives velocity and 
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temperature values which are hot too different from those observed at non-quiet 
times . This leads us to the belief that an accurate calculation could account for 
the observed values of V and T by the use of a 2 -fluid model with heating out to 
some characteristic distance R x and no heating in the interplanetary medium. 
Some preliminary calculations by Hartle and Barnes (1969) support this 
hypothesis. 

Given the T-V relation, it is possible to identify regions having anomalously 
high temperatures. We have found that these regions usually occur during large 
positive gradients in bulk speed. This is consistent with heating due to the 
collision of a fast stream with slower plasma. Heating does not occur at most 
negative bulk speed gradients, which suggests that the Kelvin-Helmholtz 
instability is not an important heating process. Heating by the action of hydro- 
magnetic shocks is not an important heating process due to the relative rarity 
of these events at 1 AU. 

V. HELIUM OBSERVATIONS 

The first long term observations of the solar wind, made by Neugebauer and 
Snyder (1966) on Mariner 2, showed that a second peak in the energy per unit 
charge spectrum was often present. This occurred at twice the energy per unit 
charge that characterized the H + peak, and was ascribed to H e ++ , an identifica- 
tion which has been confirmed by the present experiment. Thus, the solar wind 
should be considered 1 a 3-component fluid - protons, electrons and H e ++ . 

Explorer 34 observations (Ogilvie and Wilkerson, 1969) show that the bulk 
speeds of H + and H e ++ are equal on average. This is illustrated in Figure 3, 
which shows the independently determined hydrogen and helium velocities for a 
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period of 50 minutes during which they increased by 8%. The accuracy of 
determination of the helium velocity, is of course, less than that of the hydrogen 
velocity, but the two are equal to within the combined errors. 

The relative abundance of H e ++ is of basic importance in the study of the 

fluid properties of the solar wind and of basic astrophysical significance as well, 

and has been studied by several groups (Robbins et al., 1969, Neugebauer and 

++ + 

Snyder, 1966). They can accurately determine H e /H only when the tempera- 
ture is low, since they use only an energy per charge analysis. The Explorer 34 
experiment separates H e ++ and by both a mass per charge and energy per 
charge analysis and can in principle obtain H e ++ /H + for any temperature. 
However, due to an instrumental background it was often not possible to determine 
a value for the abundance n a / rip , when this quantity was low simultaneously with 
the proton density, n p , and temperature. Yet there were times of high density and 
temperature when H e densities could be measured, so observations which 
confirm and extend earlier results have been obtained. 

In Figure 4 we show the frequency distribution of all observations of n a /n p 
obtained during the life of the experiment. The number of measurements 
giving a relative abundance of zero (that is , no helium detected) is not shown and 
will be discussed below. This distribution is skewed and is not a normal 
distribution, as it would be if there were a unique constant relative abundance 
and the variability were the result of measurement errors. The presence of 
the background, although well determined and very constant, is one obvious 
cause for the skewing observed in Figure 4. We can conclude from this dis- 
tribution that the average value of n a /n p is approximately 0.05 ± 0.015, and that 
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there is considerable variability. The average value is, if anything, too high due 
to the exclusion of periods of time when the relative abundance of helium is 
especially low. The variability is a real physical effect, associated with solar 
conditions, and has time scale of order hours or tens of hours. There exist 
many continuous subsets of the data, lasting between a half day and two days 
and containing or order 500 consecutive readings, whose distributions are quite 
distinct from one another. The average value of n^np = (0.051) compares 

with the result of Neugebauer and Snyder, 0.046 ± .038, and of Strong, Asbridge, 
Bame and Hundhausen, 0.046. 

In Figure 5 we see two distributions of n a /n p , the solid histogram being 
for measurements made on June 26-27, 1967, and the dotted histogram for 
measurements made earlier on June 8, 1967. Although the latter histogram 
contains relatively fewer events, the gross difference in characteristics between 
it and that for 26th-27th June is apparent. The most probable value of n a /n p 
is higher and the spread also probably greater. 

The detection limit for the experiment is a complicated function of the 
plasma density, velocity and temperature. If the temperature is low, the bulk 
speed must coincide with a differential velocity channel or no counts will be 
recorded. If the temperature is high, enough counts must fall into each of the 
differential velocity channels covered by the distribution to exceed the background 
by twice the standard deviation of the background. For times when the proton 
density was high, (>10), the apparent absence of helium always coincided with 
low temperatures and bulk speeds falling appropriately between instrumental 
channels. Thus this biased sample does not provide evidence for an occasional 
complete absence of helium from the plasma. 
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Periods when the ratio n a /n p is greater than normal sometimes closely 
follow geomagnetic storms. Although it is tempting and almost certainly 
correct to assume that helium-rich plasma occurs as a result of flare acrivity 
on the sun, often as part of the driver of a shock, the data are consistent with a 
random association. In view of this and of known difficulties of associating storms 
and flares, we must regard this case as not entirely proven. 

The ratio of the helium temperature, T tt , to that of the protons, T p , is of 
particular interest from the point of view of understanding heating processes in 
the solar wind. If T a = T p , the two species are in thermodynamic equilibrium, 
whereas if T a = 4T p they may be characterized by the same velocity distribution 
function. As pointed out by Jokipii and Davis (1968), a process which causes a 
change in the velocity distribution function depending only upon the particle 
velocity will give ions the same velocity distribution, and temperatures propor- 
tional to their masses. Thus, observations showing that T a /T p ~4 are consistent 
with heating by hydro magnetic waves , while equal temperatures would be con- 
sistent with collisional heating leading to thermal equilibrium. Such observations 
cannot tell us, however, where the heating took place, since adiabatic expansion 
conserves the ratio. 

There have been many periods of time when more than 2 energy channels 
for the helium ions register counting rates above the background. It is then 
possible to determine independent values of the temperatures and other fluid 
parameters for each species during the same time period, by the curve fitting 
method normally used for the protons. Examples of the results of this procedure 
are shown in Table III. The periods of time indicated in the second column are 
those for which conditions remained constant; the fluid parameters set out are 
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averages of those determined every three minutes during these periods. The 
proton density n p> set out in the third column, is higher than average; this 
reflects the condition for good helium observations to be made by the p'es.-,nv; 
experiment. The two independent determinations of the bulk velocity of the 
plasma, in the last two columns, show remarkably good agreement. The ratio 
T a /T p is found to be always greater than one, and usually less than 4, suggesting 
that some form of hydromagnetic wave heating mechanism was operative. The 
helium observations thus agree well with the hypothesis of heating close to the 
sun, but cannot provide a decisive test of it. 

VI. FLOW BEHIND SHOCKS 

Shocks occurred on 11 August, 25 June and 26 June, 1967 when both 
Explorer 34 and Mariner 5 were measuring solar wind parameters at different 
locations. Examination of the fluid parameters behind the shocks showed a high 
degree of correlation between the measurements made at both separations; thus 
one may with confidence talk of the flow regime behind the shock fronts. In 
collaboration with Dr. A. Lazarus, we have examined the observed flows in 
terms of existing models: 

1) . An instantaneous coronal heat source leading to a spherical shock. 

This situation has been treated by Parker (1963), Rogers (1957), Dryer and 
Jones (1968) and Hundhausen and Gentry (1969). 

2) . An enhanced coronal heat source present [or a finite time leading to a 
shock driven by a spherically symmetric plasma. This has also been treated by 
Hundhausen and Gentry (1969), and Parker (1963). 
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3) . A shock associated with a narrow jet of plasma ejected from the site 
of a flare, as suggested by Gold (1955) and by Akasofu and Yoshida (1967). 

4) . A quasi- stationary shock aligned along the spiral direction, produced 
by a long lived, co-rotating, high speed stream. 

Hourly averages of the solar wind parameters observed on Mariner 5 were 
used to predict the parameters seen at earth by making the following adjustments: 

1) the number density was reduced by assuming an inverse square dependence 
on distance from the sun. 2) the convection speed, V, was assumed to remain 
constant, and the arrival time was obtained by dividing the difference in distance 
from the sun by that speed. These adjustments were applied to each of the hour 
averages. Note that we assumed an advancing "front” of plasma perpendicular 
to the sun-earth line. During the event of August 11, the spacecraft was essentially 
on the sun-earth line and 1.6 x 10 7 Km from the earth. During the events of 
25 and 26 June, Figure 6, the distance separating the two spacecraft was less 
than ,1.7 X 10 6 Km. 

Consider the two shocks on 25 and 26 June. Figure 6 shows that the flow 
behind them is very complicated, possibly due to the short elapsed time between 
them. During the quiet periods immediately before the shock on 25 June, and 
around 1200 hours on 26 June agreement between the two plasma instruments 
was very good. The velocities and most probable thermal speeds agreed well 
at other times, and the disagreement between the density observations for the 
period 1800 on the 25 June to 0600 on 26 June is taken to indicate that the scale 
of density fluctuation was as small as 0.01 AU. 

The bulk speed increases monotomeally during the period of observation, 
and the mean thermal speed, equivalent to the temperature, shows three peaks, 
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one immediate!}' after each shock and one between them. The increasing 
velocity and temperature behind the shocks is inconsistent with their being blast 
waves, behind which the models predict monotonically decreasing density, velocity 
and temperature. This indicates that the shocks are probably driven, but the 
driving mechanism cannot be determined by the present experiment. 

Having established the intercalibration of the two instruments, we now turn 
to considerations of the Aug. 11 event. Both spacecraft were situated on the 
earth-sun line, and separated by a distance of 1.6 x 10 7 Km. Figure 7 shows 
hour averages of the bulk speed, density and most probable thermal speed 
measured by the MIT instrument on Mariner 5, adjusted for the expected 
inverse square density-distance dependence and for distance. These are com- 
pared with the equivalent quantities observed at the earth by the Explorer 34 
instrument. The agreement of all three quantities is striking, except for the 
period between 0600 and 1800 on August 11. 

Let us now consider the bulk speed results at the top of Figure 7. The shock 
arrived at Explorer 34 at 0555 UT on 11 August, and the bulk speed increased 
monotonically until about 1600 UT. The drop at about 0800 is not real, being 
due to the deflection of the plasma out of the cone of acceptance of the instrument. 
The shock occurred at Mariner 5 during a data transmission gap between 0000 and 
0600, and the bulk speed reached its maximum at about 1300. Since the time delay 
due to the 1.6 x 10 7 Km separating the two spacecraft has been corrected for, the 
remaining time difference, seen mostly clearly in the bulk speed curve of 
Figure 7, must be due to the velocity (U-V) of the shock relative to the plasma. 
Using the times of observation at the two spacecraft we find U > 749 Km sec ” K 
This is a reasonable value for such a shock speed (Ogilvie and Burlaga, 1969) 
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although higher than usual. Th‘e very good agreement shown between the 
observations of all three quantities indicates that both spacecraft were studying 
the same material, which was moving with little distortion over the length 
scale of 1.6 x 10 7 Km, or ~ 0.1 AU. 

The data are consistent with the material seen before the velocity maximum 
being ambient gas heated by the passage of the shock and piled up in front of the 
driving gas. This period of time has been identified as a region of local heating 
(Burlaga and Ogilvie, 1969b). The shock is then to be interpreted as M standing 
off" the advancing high speed stream. The observations qualitatively resemble 
the results Hundhausen and Gentry (1969) for driven shocks. If we identify the 
stand-off distance with the product of the bulk speed and the time between the 
shock and the velocity maximum, as shown in Figure 2, we see that this distance 
is of order 0.15 AU and that it increased by about 25% in going from the position 
of Mariner to i AU. This value is somewhat smaller than that predicted by 
Hundhausen and Gentry for a spherically symmetric shock. A small stand-off 
distance and a weak shock can be simultaneously achieved by postulating that 
the advancing driver gas is not spherically symmetrical. Using the earth’s bow 
shock as an analogy, we find from the work of Spreiter and Jones (1963) that the 
radius of curvature of the front must be about 0.6 AU, This is consistent with 
observations of Akasofu and Yoshida (1967) and Hirschberg (1968) and also with 
shock observations by Taylor (1968), who noted departures from spherical 
symmetry in shock fronts. 

We thus find that although the main features of the 11 August event can be 
described well in terms of a radially propagating driver shock as calculated by 
Hundhausen and Gentry, some quantitative details favor a tongue of plasma with 
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a radius of curvature of about 0.6 AU similar to the Gold model* These 
observations are of interest because by the use of two satellites we can definitely 
rule out a co-rotating shock in this case. 

CONCLUSION 

We have interpreted approximately 3000 hours of satellite observations 
with reference to a hydromagnetic fluid model of the solar wind. We find that 
this model describes the main features of the observations and that many 
hydromagnetic phenomena are exhibited, including shocks, heating due to 
waves, and discontinuities, and also that the solar wind has a variable chemical 
constitution. These observations provide obvious pointers for the direction of 
future research into a rich variety of phenomena. 
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si’s and ssc’s Reported by 10 or More Stations, June-December 1967 
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♦Identified as a hydro magnetic shock by Ogilvie and Burlaga. (1969). 


Table II 


Ri(X 10 6 Km) 

V (Km Sec" 1 ) 

T (°K) 

T (observed, °K) 
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6 x 10 3 

1.2 x 104 
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EMPERRTURE 



Figure 1. A plot of hourly average values of the square root of the temperature 

as a function of bulk speed. 
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U (km/sec) 

Figure 2. Values of T /2 , computed from three hour averages, plotted as a function of bulk speed V for intervals of 50 km sec' 1 . The 
open dots are Explorer 34 results, and the solid dots are Parker’s solution to Bernoullis equation for an isothermal corona, TaV 2 
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is histogram shows the results of the measurements of n m /n p for the period of operation of Explorer 34. 




Figure 5. Two histograms showing observations of n a /n at different times. 
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Figure 6. Comparison between Mariner 5 and Explorer 34 Observations on June 15 and 26, 1 967. 
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